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Abstract 



> 

■ The exact nonequilibrium time evolution of the momentum distribution 
O " 

■ for a finite many particle system in one dimension with a linear energy dis- 

o 

£ — , ■ persion coupled to optical phonons is presented. For distinguishable particles 

the influence function of the phonon bath can be evaluated also for a finite 
£3 ' particle density in the thermodynamic limit. In the case of fermions the 

exact fulfillment of the Pauli principle involves a sum over permutations of 

o ' 

O ■ the electrons and the numerical evaluation is restricted to a finite number 



of electrons. In the dynamics the antisymmetry of the wavefunction shows 



5_i 

up in the obvious Pauli blocking of momentum states as well as more subtle 

interference effects. The model shows the expected physical features known 

from approximate treatments of more realistic models for the relaxation in 

the energy regime far from the bottom of the conduction band and provides 

an excellent testing ground for quantum kinetic equations. 
PACS numbers: 71.38.+i, 78.47.+p, 72.10.Bg 
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I. INTRODUCTION 



The short time evolution of the relaxation of hot electrons due to the emission of opti- 
cal phonons is not properly described using a Boltzmann equation, which assumes energy 
conserving independent scattering events, but requires the use of quantum kinetic equations 
|T|-^]. Uncontrolled approximations involved in the derivation of these equations can lead 
to unphysical results such as negative probabilities for the momentum distribution. For 
the testing of new approximations it is very useful to find exactly solvable models, which 
show the main expected properties. We recently proposed two one-dimensional models for 
the relaxation of a single electron in which the energy dispersion of the excited electron is 
linearized around the initial energy high in the conduction band. In the first model [g] an 
infinite Fermi sea, as in the Tomonaga-Luttinger-model || is included. This can be viewed 
as a crude model for the relaxation in a doped semiconductor , while the second model 
without the Fermi sea is more relevant for undoped semiconductors. Exact results for 
the momentum distribution were presented and compared with various improved quantum 
kinetic equations. 

In a typical experimental situation a finite density of electrons is optically excited into 
the valence band. The nonequilibrium momentum distribution can be inferred from energy 
resolved measurements of transmission changes AT/T with delay times in the femtosecond 
range jlO |. For their theoretical discussion of the results for GaAs these authors have com- 
pared approximate calculations in more realistic band structures with our second model @ 
and "find the ID linear model to be well justified as long as the carriers do not reach the bot- 
tom of the band" . For a more detailed study of the relaxation process it is necessary to treat 
the two effects properly which are not correctly described if the results of the single electron 
model are simply superposed according to the density of the excited electrons. A phonon 
emitted by one electron can be reabsorbed by another one leading to energy gain satellites 
in the momentum distribution. In the relaxation of many electrons the Pauli principle has 
to be taken into account. In a Boltzmann equation or typical quantum kinetic equations the 
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Pauli blocking of states is only treated in an average way and no exact description has been 
provided previously. In this paper we generalize our second model [[!]] to the many electron 
case and again obtain the exact solution for the nonequilibrium momentum distribution. 

In Sec. I| we present the model and its solution for two cases. We first assume the 
"electrons" to be distinguishable. Under this assumption the results for finite systems with 
a finite number of particles can be extended to the thermodynamic limit with a finite density 
of particles. The results show the importance of the gain processes mentioned earlier. If 
the fermionic character of the "electrons" is properly taken into account, the result for 
the momentum distribution contains a sum over permutations of iV — 1 electrons and a 
numerical evaluation of the exact result is only possible for finite systems. For not too large 
electron-phonon coupling we also present an excellent approximation to the exact results 
which can be evaluated for larger systems. The model allows a detailed study in which 
parameter regimes a proper treatment of the Pauli blocking is of importance. Interference 
effects which are not correctly described by the usual quantum kinetic equations are shown 
to be especially important in the gain region. Some analytical and numerical details of the 
solution as well as a comparison with time dependent perturbation theory are presented in 
appendices. 



II. THE MODEL AND ITS SOLUTION 

In generalization of the polaron model presented in we study a one-dimensional model 
of the relaxation of iV "hot" particles due to the interaction with optical phonons. In order 
to obtain the exact solution of the time dependent Schrodinger equation we assume a linear 
energy dispersion for the particles. The Hamiltonian for a finite system of length L with 
periodic boundary conditions reads 

H = H e + Hp + H ep = H + H ep (1) 

with 
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N 



i=1 1/2 9 ( 2 ) 

where t> is the velocity of the particles and b q the creation operator of a phonon with 
frequency All the numerical results presented in Sees. [TT1] and [TV] are for optical phonons 
with a constant phonon frequency u q = ujq. The "electron-phonon" coupling strength given 
by g(q) is specified later and the operator p q denotes the Fourier transform of the particle 
density 

N 

ft = E^- (3) 

1=1 

We will discuss the case of distinguishable particles as well as the model with N spinless 
fermions. Due to the linear dispersion the model has no ground state. In the fermionic case 
this could be altered by adding an infinite Dirac sea. This is not considered in the following 
and we study the relaxation of an initial state with initial momenta k®, . . . , k% and the 
phonons in their unperturbed ground state or in thermal equilibrium at a finite temperature. 
The exact description of the time dependence is possible due to the commutation relation 
[£i, H] = ivH, which leads to 

p q (t) = e iHt p q e- iHt = e-^p, (4) 
independent of the coupling strength g(q). The Heisenberg equation for b q (t) reads 

ih q {t) = u q b q {t) + C^f] ' g{q) Pq {t). (5) 



dt BV ' q qK ' \L 
The solution of this inhomogeneous equation follows with 



b q (t) 



e 



-iiOqt 



(6) 



Therefore the change of the expectation value of the number of phonons N q = b q b q due to 
the presence of the particles is determined by the expectation value of p q p q in the initial 

state 
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6{N q ) {t) = ^\f q {t)\ 2 (p\p q 



(7) 



with 



(8) 



This will be discussed in Sees. |TTT| and [TV]. We first present the method to calculate expecta- 
tion values of an operator A acting in the Hilbert space of the particles, i. e. A = \k)^ ^ (k\, 
where A = \k)^ is the momentum state in the Hilbert space of the particle i. For the case of 
indistinguishable particles we will consider one-particle operators A = £)i=i A^. We present 
the result for expectation values (A) (t) for a statistical operator of the initial state which 
factorizes into a particle and a phonon part, i. e. Po = p p ■ p p h- We work in the interaction 
representation and denote the time evolution with the unperturbed Hamiltonian by Ao(t). 
Due to the simple form of Pq : D(f) — e~ iqvt p q it is favorable to perform the trace in the particle 
Hilbert space using the position states |x) = (g) • • • ® \xn) n . This yields 



(A) (t) = / dx'dx (x' \p p \ x) Tr 



ph 



Pph <x 



&{t)A D (t)U{t) 



(9) 



where U(t) = e tH ot e -tHt j g ^ e time evolution operator in the interaction representation. The 
proper projection on the totally antisymmetric states in the case of fermions is taken care 
of by the statistical operator p p . In the following we assume p p to describe a translational 
invariant initial state, i. e. p p is diagonal in the momentum representation. The simple form 
of p q ,D{t) implies that the only electronic operators in H ep ^{t) are the commuting position 
operators Therefore U(t) |x') = U x r(t) |x'), where U x r(t) is an operator in the phonon 
Hilbert space only. This allows us to write Eq.(0) in the form 



(A) (t) = J rfx'rfx (x'|p p |x) (x\A D {t)\x!)F(x.,x/,t), 



(10) 



where the "influence function" -F(x, x',t) is given by 



F(x,x',t) = Tr ph \p ph Ui(t)U x ,(t) 



(11) 



The differential equation for [7 x (i) reads 



dt 

with 



du x {t) = (Cjf) + cl (t)) U x (t) (12) 



CxW = (^J E9(9)^ r " ? V ? W6 9 , (13) 

where p 9 (x) is the eigenvalue of the operator p q in the state |x). Using the Baker-Hausdorf 
formula the solution of equation (|P2"D is straightforward [0] 

l/ x (*) = exp | -i J Cl(t')dt'} exp |-i J* C x (t')dt'}a x {t) , (14) 

where a x (i) is given by the c-number 

= exp {-|ET ? W Ip,(x)| 2 } (15) 



with 



%(t) = . g2{q) , 2 (l - i K - f - e-^^)*) . (16) 
[u q -vq) K 



The trace over the Hilbert space in Eq.(|TTJ) can be performed if the phonons are initially 
assumed to be in thermal equilibrium at temperature T. Then we can use (e A e B 

1 / j2_|_9 d d_j_ d2\ 

e 2 \ /, valid for canonical expectation values for free bosons and operators A and 

B linear in the boson creation and annihilation operators M. Using Eqs.(|T3|) and ( |T4| ) one 
obtains 

F(x,x',i) =exp^J2\fM 2 [ff q (x)p q W-n B {u q ) |p,(x) - p 9 (x')| 2 ] }<(t)a x ,(t), (17) 

where nB{ui q ) = 1/ (exp(u g / ksTph) — 1) is the Bose function. With this analytical expression 
for the influence function the evaluation of expectation values of particle operators is reduced 
to integrations in Eq . (|lOl) . 

We now specialize the operator A to the momentum distribution and begin with dis- 
tinguishable particles. For the operator rv^ = \k}^ (i) (k\ the integral in (|10|) simplifies 
as 
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(x|4 1] |x'> = ^^-^^(x-x) 



(18) 



contains a delta function in the variables x = (x%, . . . ,njy) and x'. Therefore the influence 
function F is needed only for x' = x. In the variables has 

N N 



j=2 i,j=2 



(19) 



and corresponding expressions for |p g (a;i,x)| 2 and \p q (x' l ,~x}) ] \' . Due to the identity jg(t) + 
j q (t) = \f q (t)\ 2 the terms with the double sum involving two "other" particles cancel. This 
leads to 

F{ Xl , x; x[, x, t) = F 1 {x u t) f[ exp j ^ £ [ 7 *(t) (e^ 1 - l) e"** - lg (t) (e 4 ** - l) e 1 ^] J 

(20) 



where 



Fi(x, t) = exp |^ 2 l/«(*)| 2 [( e ^ - X ) C 1 + + ( e ~^ - l) ^K)] } 



(21) 



is the influence function of a single particle [|7| . 

Now the integration in Eq . (|i0|) can be partially carried out using Eqs. (|i"8"D and (|20|) and 
one obtains 



n «\ (t) = [ dxe ikx p^ 1 \x)F 1 (x,t) (G (x,t)) 



N-l 



(22) 



with 



1 

I 



G (s, f) = i j dy exp { ^ £ [ 7 *(t) (>* - l) e'^ - lq (t) (e 



-iqx 



- 1 e tqy 



(23) 



and p^ (x) is given by the real space matrix elements of the reduced statistical operator 
for particle one for which the trace (tr r ) over the other particles has been performed 



P {1 \x l -x' x ) 



N-l 



(xi,x| p p K,x) = (xi\tr'(p p ) \x[) ■ 



(24) 



For N = 1 and p^> = \k®) , i. e. p^(x) = \e lk i x Eq.(^) reduces to the polaron result 
presented in |7J . The result for (njj. \ (t) will be analyzed in the next section. 



We now turn to the case of N fermions. In first quantization the momentum occupation 
number operator n k reads in the Hilbert space where (a) denotes the subspace of 

antisymmetrized states 

njfe = jj-j— \k,k 2 , . . . ,k N ) aa (k N , . . . ,k 2 ,k\ (25) 

The factor 1/(N — 1)! is necessary because the N — 1 summations are unrestricted. To 
calculate (h}~) (t) we can proceed as in Eqs.(|9]) and fllCf ) using product states |x) which are 
not antisymmetrized. In the following we describe the time evolution of an initial state 
| ko) . The result for an arbitrary translational invariant statistical operator is obtained by 
averaging the result with the normalized probability density p(k ) describing p p . The initial 
state \ki, . . . , k° N ) a and the projectors in Eq. (|25|) take care of the Pauli principle. Using 
(x|k) a = a (x|k) we obtain 

where the parenthesis contains all (iV!) 2 terms which result from antisymmetrizing the first 
term with respect to the x and the x' variables. As in Eq. ([T0|) this expression is multiplied 
by the function /(x, x',t) = (k°|x) aa (x'|k°) F(x, x',t), which is also antisymmetric in both 
the x and x' variables, the x— and x'— integrations yield (N\) 2 identical terms, i.e. one can 
just take the first term on the rhs of Eq.(|26|) and multiply it by (N\) 2 

(n k ) = — J da; 1 dx , 1 tface tt(xi - x i ) o (k |x 1 ,x) (a;;, x|k°) a F(xi,x; x'^x,*). (27) 

Note that again only the special form of the arguments of the influence function as in 
Eq.flTTD appear. Due to the translational invariance discussed there, the x' x — integration can 
be carried out by setting x[ = in the remaining integrand and dropping the factor 1/L. 
The factors in the integrand which take care of the Pauli principle can be written as 

i N 

a (k K,x) (0,x|k ) a = — £ e ifc ^ a (k | C lo|x)(x|c fc o|k ) a (28) 

where the c\,(ck) are the usual creation (annihilation) operators and the factor 1/N appears 
because we have switched from antisymmetrized iV— particles states to antisymmetrized 
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C '^(x-; 



± 



(26) 



(iV— 1)— particle states. In Eq. fl2"7|) the rhs of Eq. fl28|) is integrated with a function symmetric 
in the x variables. In writing out the factors a (k°| c k o |x) and (x| c k o |k°) a as antisymmetrical 
sums of (iV — 1)! terms each term of the second sum yields the same contribution. Therefore 
for the use of the integration with a symmetric function in x we can write 



a (ko|xi,x) (0,x|k ) c 



NL N 



N 

E 



J(k-k°)x 



1 {-iy +j+PN - i e 



(29) 



»J=1 Pn-i 

Here we have introduced the notation (1 = 2,..., N) 



j = 



k? (I > J) 

i.e. the upper index indicates the missing initial momentum. If one inserts 
of Eq. (|27D one obtains 



(30) 



on the rhs 



N i 

(n k ) (t) = Et / dxe^^Fi^^GSli^t) 

8=1 



(31) 



with 



N 



N 



(32) 



1=2 



and 



G(k,x, t) = jj dye^exp £ [t?(*) - l) - 7 ,(t) (« 



-1 e 



(33) 



This is a central result of our paper. The result for (n k ) (t) in Eq.([JT]) should be compared 
with the result for (n k ^ (t) in Eq. (|22|) summed over all iV particles. If all the initial momenta 
are different also for the distinguishable particles the total momentum distribution is given by 
Eq.(|3~TD with G^^x, t) replaced by (Gq(x, t))^ -1 . For fermions the sum over permutations 
in Eq.(^) is necessary in order to obey the Pauli principle exactly. While for distinguishable 
particles it is possible to obtain results for a finite density of particles in the thermodynamic 
limit as discussed in Sec.|T|, the more complicated expression for G^ ) _ 1 (x, t) can only be 
evaluated for a finite number of particles. The results are presented in Sec. |V|. 
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III. RESULTS FOR DISTINGUISHABLE PARTICLES 



As we are mainly interested in modelling aspects of hot electron relaxation in semicon- 
ductors, the fermion result Eq.(|3T|) for the time dependence of the momentum occupation 
numbers is considered more important than the expression Eq.(p2|) for distinguishable par- 
ticles. Due to the relative simplicity of the latter result compared to the fermionic one, it is 
useful to discuss Eq.(^) in order to obtain a first understanding of the difference between 
the cases N = 1 and N > 1. 

The modification of the energy transfer to the phonons is trivial, as ^k p\p q k ^> = N 
for q 7^ , i. e. the N = 1 polaron result is just multiplied by N in Eq. ([?]). In contrast 
to this simple expression, the occupancy (^i ) (t) for N > 1 differs qualitatively from the 
N = 1 result. With the phonons initially in their ground state no energy gain is possible for 
one particle, if in our model the g-sum in H in Eq. ([I]) is restricted to positive values. If the 
sum also includes negative g-values a small weight in the gain region occurs for short times 
but no resonant gain peak appears. In the following we therefore use two special choices 
for the coupling function of the "electron" -phonon interaction, g^\q) = gQ(q)Q(q c — q), 
where q c is a cutoff and g 2 (q) = g for all q G (—00,00). As discussed in appendix |A| 
the latter choice allows the analytical calculation of the integrals in the influence function 
-F(x, x',t) in Eq. (|T7|) in the thermodynamic limit. We start our discussion with the special 
case p^ = tr'p p = \k\) (k®\ , i.e. particle one is initially in the state \k®). 

A straightforward perturbation expansion of Gq(x) in Eq.(|23D in powers of 7 ~ g 2 shows 
that energy gain is possible for N > 1 in contrast to the N = 1 result for g^(q) : 

Go(^,t) = l- (^) 2 E{[l7,W| 2 + ^( 79 (t) 7 _,(t))] (2-e^-e-^)} + 0( 7 3 ) (34) 

For the model with g^(q) the contribution proportional to e~ iqx leads to a nonzero 
weight in (n^^ (t) for k > k®, i.e. energy gain. Only the term proportional to |7g(£)| 2 is 
nonzero and the probability for the gain process agrees with the result from squaring the 
amplitude in second order time dependent perturbation theory for the process in which 



10 



particle 2 emits a phonon, which is reabsorbed by particle 1 (see appendix As |7g(£)| 2 is 
of the order (1/L)° the probability of such a gain process goes to zero as 1/L for L — > oo, 
if only two (or a finite number of) particles are present. It is therefore more interesting to 
study the thermodynamic limit for a finite (other) particle density n p = (N — 1)/L. In the 
limit L —>■ oo the momentum distribution goes over to a function defined on 1R. In order 
to normalize this function to one, (n^ \ (t) has to be multiplied by the density of k-points, 



i. e. we consider the function n(k,t) = lini^oo ( n fc ) As the function Go(x,t) in 



Eq.(p3|) has the form 



1 



G (x,t) = 1 + — _ 1 ^ p c L (x,t) 



(35) 



where ci{x,t) goes over to a well defined limiting function c(x,t) for L — > oo, the fac- 
tor (Gq(x, t)) N ~ l in Eq.(^) has the iV — > oo limit exp (n p c(x, t)). This yields with the 
normalized probability p(/ci) for the initial momentum the exact finite density momentum 
distribution in the thermodynamic limit 



n(k, t) 



dk\p{k\ 



with 



b(x,t) 



dq\f q (t)f 



jqx 



OO Ar*> 

e i(k—ki)x b(x,t)+n v c(x,t) 

-oo 2tt 



l) (l+n B K)) + (e-^-l) n B (u q ) 



(36) 



(37) 



and 



c(x, t) 



dy 



exp 



d Q 7o(*) ( e 



1 e 



7g 



-tqx 



1 e 



"/;/ 



(38) 



The function c(x, t) cannot be determined completely analytically even for the coupling 
function g { - 2 \q) = g, as discussed in appendix [A]. It is therefore useful also to consider c(x, t) 
in leading order perturbation theory in 7. Using Eq.(|34]) or (^) one obtains 



c 2 (x,t) = -2tt / dq \ lq (t)\ z + Re( lq (t) 1 _ q (t)) 



2-e 



■tqx 



- e~ iqx ) 



(39) 



In this approximation the result for n^(k,t) has the same form as for a single particle and 
the phonons at finite temperature f7j, but with the Bose function in Eq.([37|) replaced by the 
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sum of riB{uJ q ) and the time dependent quantity 2im p |7 9 (t)| 2 + Re ( r y q (t) r y^ q (t)) / \f q (t)\ 2 . 
In order to estimate in which parameter range Eq.(^) provides a good approximation it is 
useful to introduce the dimensionless coupling constant a = 2ng 2 / (vu ) and to multiply n q 
by the relevant length v/ujq = 1/qs in our model. If in addition we introduce dimensionless 
time and space variables r = u^t and u = qsx, the finite density contribution n p c(x, t) in 
the exponent in Eq.(^) can be written as (n p /qB)c(u, r) where c is a dimensionless function. 
The small density regime is given by h p = n p /qB <C 1- For finite systems this corresponds 
to (N — l)/(27rns) "C 1, where ns is the number of momentum states in an arbitrary 
momentum interval (q, q + The use of c^(x,t) in Eq . (|36D is certainly allowed in the 
low-density and small a regime. 

For the special case of a sharp initial momentum p(ki) = 6(ki — k\) the momentum 
distribution n(k,t) in Eq. (|36|) consists of a delta peak at k = k® with a weight Pfco(t) = 
exp (6(00, t) + n p c(oo, t)) and a continuous part which for finite times is nonzero for k\ also. 
In Fig.|I| we compare the time evolution of the weight p k o (t) for the two different values of 
the dimensionless density h p = 10~ 3 and n p = 10~ 2 with the N — 1 polaron result. The 
dimensionless coupling constant has the small value a = 27rl0~ 3 . The result for the coupling 
functions and do not differ on the scale of the figure. Obviously the presence of the 
other particles accelerates the relaxation process. In order to obtain a simple quantitative 
estimate of this effect we compare the exact result with the approximation to replace c(oo, t) 
by 02(00, t), where 02(00, t) is obtained by dropping the a;— dependent part on the rhs of 
Eq . (|39D . In this approximation the exponential decay of the polaron case is replaced by the 
faster decay approximately given by Pfco(r) »s exp(— ar — | a 2 h p r 3 ), where we have used 
Eq.( |A10| ). For dimensionless times r larger than («n p ) -1 / 2 the deviation from the polaron 
result becomes important. 

In Fig.0 the first two curves represent the continuous part of the momentum distribution 
as a function of the dimensionless momentum variable k = k/qs for k® = and the weight 
of the delta peak is given in the figure captions. We compare the results for the coupling 
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functions and for the same values of the coupling a as in Fig.|T] and n p = 
at the time r = 50 where the influence of the other particles in the energy loss regime 
at negative momenta is still rather small. The momentum distribution shows loss features 
around k — — 1, —2 and —3 corresponding to the emission of one, two and three phonons. 
The short time energy uncertainty well known from standard derivations of Fermi's golden 
rule is clearly visible, especially at the tails of the one phonon loss peak. For longer times 
these oscillations weaken and the peaks become narrower at the resonant positions. A 
rather small energy gain peak at k — 1 is also visible in FigfJ. A qualitative difference 
of the coupling functions and g^ shows up at positive momenta. While the weight is 
strictly zero for h v = in this range for g^(q), the coupling g <y2 \q) shows a weak oscillatory 
weight there as in the polaron case which survives in Fig.^|. Within the plotting accuracy 
the results in this figure agree with the exact results, when c(x, t) is approximated by C2(x, t) 
(Eq.(|39|)). For much larger times many sharp phonon replicas develop and the oscillations 
in the tails become weak. The difference between the g^ and g^ results become very small 
and the replacement of c(x, t) by C2{x, t) is quantitatively no longer sufficient especially in 
the gain region. The oscillatory behaviour of both systems is smeared out and the peaks are 
broadened if a Gaussian probability distribution p(ki) is used in Eq.(|36|). 

In Fig.|3| the momentum distribution functions for the thermodynamic limit with different 
particle densities and for a finite system are compared at the time u>ot = 30. Our method to 
perform the integration in Eq. for finite systems is described in appendix [TJ. Only the 
zero momentum state was initially excited but due to the larger coupling a = 27rl0 -2 the 
curves are smoother and the relaxation faster. Also, for the first two systems with the larger 
particle density the energy uncertainty oscillations of the first resonant state are smeared 
out. They have the same dimensionless particle density, i.e. for the thermodynamic limit it 
is rip/qB = 0.2/(27r) and the finite systems n p /qB = (N — l)/(27m/3) where qs = ^j; n B- As 
can be seen the result for the finite system with nine particles is already very close to the 
result for the thermodynamic limit. At earlier times the two curves agree very well but then 
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the differences appear and grow with time advancing. For a smaller coupling this behaviour 
is slowed down, e. g. for a = 27rl0 -3 the agreement is still excellent at time u^t = 50. 

Comparing the first and the last system one can see that the gain satellites at positive 
momentum and also at momentum less than — are much more pronounced for the larger 
density. Generally, a larger density smoothes the curves and accelerates the spreading of the 
peaks. 

IV. RESULTS FOR FERMIONS 

The role of the Pauli principle for more than one fermion is quite easily seen for the 
number of phonons created Eq.(^), while the momentum distribution Eq. (|3l|) requires a 
much larger numerical effort than in the case of distinguishable particles. 

In order to obtain S (N q ) (t) we have to evaluate the expectation value of p q p q in the 
fermionic initial state, which was taken as the Slater determinant |k°) a in Sec. f|. The 
expectation value can either be calculated using Wick's theorem using the method of second 
quantization or using first quantization as in Sec. [TT| For q ^ one obtains 

{p\p*) t=Q = E M (t = 0) (1 - (n k ) (t = 0)) 

k 

N 

= N- E S k°,k° +q , (40) 

where the first equality holds for an arbitrary initial Slater determinant. The second equality 
shows that the Pauli blocking on the energy transfer to the phonons is most prominent in the 
artificial state in which all the initial momenta are separated by qs- Then (p\ B Pq B / — 1 and 
due to the Pauli blocking the number of phonons corresponding to the energy conserving 
transitions is reduced drastically compared to the case of distinguishable particles. 

In order to calculate the momentum distribution Eq. (|3l|) the functions G^_ 1 have to be 
evaluated, which involve products of the functions G(k,x,t) for different values of the mo- 
mentum argument. We first discuss the behaviour of G(k, x, t) in leading order perturbation 
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theory in the electron-phonon coupling g. Expanding the exponential function in Eq.(|33|) 
we obtain 

G(k, x, t) = 4,o + (^) (im - 7-*(*)) (e lkx - l) 

+7j(*)7*- ff (*) (l + eikx ~ ei{k ~ q)x ~ ei9X ) / 2 

+l q {t)l-k- q {t) (l + e lkx - e^+ti* - e- iqx ) /2 

+0( 7 3 ). (41) 



In the sum over j on the rhs of Eq. fl32|) the term i = j and the terms j ^ i be- 
have differently concerning a perturbation expansion. We first discuss the diagonal term. 
In the sum over the permutations the identity presents the leading approximation given 
by (G (x,t)) N ~ 1 , i.e. this contribution to G^^x, t) is identical to the exact result for 
distinguishable particles. The contributions of the other permutations are of the order 
(7/L) m G v " 1 " m with m > 2. For i ^ j there is no term in which all momentum differences 
vanish. The leading term in j/L is given by —G(k® — k®,x,t)(G (x,t)) N ~ 2 . Just keeping 
the leading order terms for the diagonal and the non-diagonal terms provides the following 
approximation for G^-i 

G$_i « (G (x,t)f- 2 \G (x,t) - £ G(/c° - k%,x,t)\ (42) 

with an error of order (7/L) 2 . Therefore in this approximation the Pauli principle is strictly 
obeyed to order g 2 . To this order (n^) (t) in Eq.([H|) can easily be calculated using Go ~ 1 
and Eq.(|5) in Eq.® 

(t) - E (i-jE j + f \h2-k(t)\ 2 } 

+~r X] |/fc?-fe?( £ )| (^,fc? _ s k,k°) (43) 

where the second line on the rhs via the approximate G^_ x takes care of the Pauli principle. 
This is most easily seen for occupancy of an initial momentum state of momentum k® 
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n k of ] (t)*l- 2 -f £ \f q (t)\ 2 - (44) 

The electron with initial momentum kf cannot scatter into the momentum k® and none 
of the other electrons can scatter into momentum kf. 

In the model with the coupling g^ l \q) gain contributions first occur in order g 4 when 
the phonons are in their ground state initially. For k larger than the initial momenta, i.e. 
k > k®, i = 1, . . . N it is therefore illuminating to expand the exact solution to this order in g. 
One can distinguish two contributions to order g A . First the approximation given in Eq.(|42|) 
can be expanded to this order using Eq . (^4|) and Eq.([|l]). The correction to the result 
for distinguishable particles in Eq . (f42|) is obtained by using (G (x, i)) N ~ 2 G{k® — k Q -,x,t) « 



G{k® — kj,x,t) and Fi(x,t) ~ 1. Only the term proportional to e~ tqx in the second line on 
the rhs of Eq.([|T]) contributes to the gain. The x-integration in Eq.(|3TD is trivial and one 
obtains 

(nk) (2) (t) = (t^E |7 fc - fc o(t)-7 fc - fc o(t)| 2 + 0( 7 3 ). (45) 



There are two types of terms missing in Eq.(^2|) in order to make G^_ x correct to order 
g 4 . They are due to single pair permutations in Eq.(|32D 



G^f = -(G (x,t)) N ~ 3 J2 (G(kl-k m ,x,t)-G(k°-klx,t))G(kl-klx,t). (46) 

Expanding this expression to order r j 2 and adding it to the approximation Eq.(]4"2|) for 
Gpf_i one obtains for k larger than the initial momenta 

K) (2) (t) = (^f) 2 J2'h-k9(t) - lk _ k o(t)\ 2 + 0( 7 3 ) (47) 

where the prime indicates that only the pairs contribute to the sum for which k^+kj—k 
is different from the other initial momenta. This result which correctly incorporates the Pauli 
blocking can also quite easily be derived in second order time-dependent perturbation theory 
for the initial state (see appendix |Cp and it is very helpful for the discussion of the results 
presented in the next figures. 

16 



In the following we present detailed numerical results for finite systems using the exact 
expression for G^ ) _ 1 (x,t) as well as the approximation presented in Eq.(^) and compare 
the results with the corresponding ones for distinguishable particles with the same initial 
momenta. 

In Fig.^| the results for the distribution function of three different systems are compared 
at time u>ot = 50 with the coupling a = 27rl0 -3 . In the first system there are six indis- 
tinguishable particles (fermions) initially occupying adjacent momentum states from — 5Ag 
to zero. In the second system the total momentum distribution for six distinguishable par- 
ticles initially occupying the same states is presented. The third system corresponds to 
two fermions in adjacent momentum states with the same dimensionless density N/^Lq^) 
as the six fermion system. For initial momenta where all the particles can relax into the 
lower resonant momentum state by emitting a phonon with momentum qg the difference 
between the fermions and the distinguishable particles is generally small in the loss part of 
the distribution up to intermediate times as can be seen in Fig.|j. The occupation numbers 
in the region of the gain peak (k = qg), however, differ roughly by a factor of two. For 
shorter times the deviations are more drastic. In appendix y the occupation numbers in the 
gain region are calculated in second order time-dependent perturbation theory for the initial 
state for fermions and distinguishable particles and the results derived there can account 
very well for this difference. The physical content of Eq . (|C5|) (agreeing with Eq . (|47D ) and 
Eq. (|U6|) is that a particle can be scattered into a higher momentum state with an amplitude 
1q{t) by absorbing a phonon with momentum q which was emitted previously by a different 
particle. This process has an exchange process where the initial and the final state are the 
same but the role of the two particles as an emitter and an absorber of the phonon are 
exchanged. The phonons exchanged in the two processes are different but this is irrelevant 
for the calculation of the fermionic momentum distribution. For fermions the occupation 
probability for a momentum state in the gain region is the absolute square of the sum of 
the amplitudes for the two processes which have a relative minus sign. For distinguishable 
particles on the other hand the two processes can be distinguished and the probability is the 
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sum of the absolute square of the amplitudes for the individual processes. In Fig.f| all the 
initially occupied momentum states are very close and consequently all the corresponding 
exchange processes have a similar amplitude leading to the reduction of the weight in the 
gain peak for the fermions. A more detailed discussion of j q (t) is given in appendix |C| and it 
shows that the shorter the time the broader are the features of the amplitude and therefore 
the cancellation in Eq. (f47|) . 

The interference effect contained in Eq. (f47D is not properly described by the usual quan- 
tum kinetic equations in which the collision term is treated in the Born approximation 
but requires an improved treatment ||. 

The agreement of the two fermion system with the corresponding six particle system is 
very good - also for larger times. The momentum grid of the smaller system is different 
but the results for its momentum states nearly always agree with the results for the larger 
system. The approximation presented in Eq.([4^) agrees with the exact solution within the 
plotting accuracy. 

In Fig.[5] the distributions of the same six particle systems as in Fig.|] are presented at 
time u t = 50 but with the coupling strength a = 27rl0~ 2 . Consequently the relaxation 
process is more advanced and the peaks are smoother with a larger weight in the states 
between them. The approximation (Eq.([|2|)) for the six indistinguishable particles - not 
shown in the figure for its clarity - gives values slightly too large for the central excitations 
and phonon peaks smaller than — 7qB and larger than 2qs- 

In Fig.|6] the results are shown for the same systems at the same time as in Fig.|5] but 
again with the weaker coupling a = 27rl0~ 3 and different initial momenta. Here the zero 
momentum state and the next two lower states are occupied and the second group of three 
particles has momenta shifted by — Therefore in the six fermion system the three particles 
with larger momenta cannot relax by emitting a resonant phonon. The distinguishable 
particles of course are not subject to this blocking and in contrast to Fig.|] one can clearly 
see the large difference between the two systems in the negative momentum peaks. However 
the difference in the first gain peak is not as large as in Fig.|j. For an explanation Eq.(f47|) 
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is useful again. If one particle of the smaller and larger momentum group interact by 
phonon exchange the corresponding exchange process has a very different amplitude. In the 
one process a phonon with momentum qs will be absorbed and in the exchange process a 
phonon with momentum 2g#. Consequently the contributions of the two processes do not 
cancel each other as much as in Fig.^. On the other hand processes involving two particles 
from the larger momentum group only do not contribute to the resonant states of the gain 
peak at all but the exclusion of these processes cannot compensate the other effect and the 
difference between the fermions and the distinguishable particles is reduced. 



V. SUMMARY 

We have presented the exact analytical result for the time dependence of the electronic 
momentum distribution for a model with a linear electronic energy dispersion for a trans- 
lationally invariant initial state with N electrons and thermal phonons at temperature T p h. 
In order to show as clearly as possible the difference between the N — 1 "polaron case" and 
the new effects occurring for iV > 1 we have presented numerical results for T p h = only. 
While for distinguishable particles the influence function can be simply evaluated also for a 
finite density of particles in the thermodynamic limit, for fermions numerical results were 
presented only for finite particle number and system size. As shown in Fig.^for not too large 
times finite systems with different N but the same density are very similar, which indicates 
that going to the thermodynamic limit is not very important. Apart from the obvious effects 
of Pauli-blocking interference effects play an important role which are not incorporated in 
quantum kinetic equations which treat the scattering term in the Born approximation |J 
and 0. As interference effects are very important to correctly describe the gain process in 
the short time limit our exact results can serve as an excellent testing ground for improved 
quantum kinetic equations. 
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APPENDIX A: 



In this appendix we present analytical results for various integrals appearing in Eqs.([22l) 
and (p3|) for the model with a g-independent coupling function g^ 2 \q) = g and u q = tu . We 



first consider the integral 

i( X , t) =r d J 2 e * ( Ai) 

which is already needed for the N = 1 "polaron case" . Changing the integration variable it 
can be written as 

I(x, t) = -e^ x r du^P-e'^ (A2) 

V J-oo U 

with y = 2x/(vt). Using contour integration one obtains for t > 

I( x , t) = ^e iqBX (vt - \x\) e(vt -\x\), (A3) 

v 

where 9{x) is the step function. Similarly the function 

J{x, t) = f°° dq- — -2 (l + % (u - vq) t - e^"-^) e iqx (A4) 

J-oo [uj — vq) v ' 

which is needed in the calculation of G (x, t) in Eq.(|23|) can be evaluated by contour inte- 
gration, too. For t > one has 

J( x , t) = ^e iqBX (vt - x)Q(x)Q(vt - x). (A5) 
v 2 
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With these results the q— integration of the functions in the integrand of the momen- 
tum distribution Eq. fl36|) can be performed analytically for g^(q). As a function of the 
dimensionless momentum variable k = (k — k\)/qB one obtains for n(k, r) = qBn(k,t) 

rOD du 



n(k, t) 



where b and c for r > are given by 



2vr 



gifcn ^b(u,T)+h p c(u,r) 



(A6) 



6(u,r) = -ar(l + 2n B (o;o)) + a e iu (1 + n B (cu )) + e~ iu n B (uJo) (r - \u\)Q(r - \u\) (A7) 



and 



c(u, t) = dv expj — 2ia[(r — {v — u)) sin(t> — u)Q(v — u)Q{t — (v — u)) 
-(r - u) sin(w)e(w)e(r - u)]} - 1 



(A8) 



The functions obey the relations b*(u, r) = b(—u, r) and c*(w, r) = c(—u, r) which guar- 
antee that n(k,r) is real. The first term on the rhs of Eq. ([A7| ) is equal to 6(oo, r) which 
enters the weight Pfco(r) and c(oo, r) is given by 



c(oo, r) = 2 / (if [cos {2a(r — t>) sin v} — 1] 

JO 



(A9) 



This integral cannot be performed analytically. Expanding in powers of a yields 

p k o(r) =exp|-ar(l + 2n B (cj ))-a 2 n p ^r 3 -r+isin2r^ +0(/)| (A10) 

Due to the form of the functions b(u, r) and c(u, r) the continuous part of n(k, r) can be 
expressed as a finite integral over u from zero to r 



n con t(k,T) = 2Re \ j —e 



i(k-k°) u 



e b(u,r)+n p c(u,r) _ p ^ 



(All) 



APPENDIX B: 



In this appendix we shortly describe our numerical method to calculate the momentum 
distribution for finite systems. It is a straightforward generalization of the technique used 
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in [0] . As one step of the procedure can actually be simplified considerably compared to [|7| , 
we first describe the method for the simplest case N — 1 , the coupling function g^(q) and 
the phonons initially in their ground state. In this case the function Fi(x, t) in Eq.([H]) can 
be expanded in a power series in z = e l( *^ x 

oo 

F 1 (x,t)= £ F m (t)e im ^ x . (Bl) 

m=0 

Recursion relations to determine the coefficients F m (t) are discussed below. If we take 
p^'(x) = 1/L, i.e. the initial momentum of the particle to be ko = 0, the integration in 
Eq.(p2|) can be trivially performed and yields for k n = n(^) 



(i) 



1 „ oo 

>(*) = T / E F m (t)e« n+m ^dx 

L J m=0 



= F. n (t) (B2) 
for n < and (n^)(t) = for n > 0. In order to calculate the coefficients F m (t) we write 

F 1 (x,t) = e gt{z) - 9til) = Be 9t{z) (B3) 
where the function gt(z) follows from Eq . as 

oo 

g t {z) = Y J dn{t)z n (B4) 

n=0 

with a n (t) = (^)\f n (2iL)(t)\ 2 . We now expand e 9t ^ in a power series in z 

oo 

F t (z)=e 9t ^ = Y / Cm(t)z m . (B5) 

n=0 

If we use F{(z) = g' t (z)F t (z) we obtain the recursion relations 

c m {t) = —Y\ Icm-iai (t) . (B6) 

m l=l 

Using F m (t) = Bc m (t) yields the momentum distribution for the simplest case. Already 
relaxing the assumption T p h = leads to a new aspect as Fi(x,t) can now be written as 
a product of a power series in z and a power series in 1/z. Multiplying these series and 
performing the integration as in Eq.( |B2[ ) leads to an expression for (rik n ){t) which involves 
a summation over a product of the expansion coefficients . 
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In order to calculate the momentum distribution for finite systems and iV > 1 we also 
express the functions G(n(^-),x,t) as Laurent series in z. As t) is a sum over 

products of these functions it also can be obtained as Laurent series in z such that the final 
integration in Eq. (plf) can be performed trivially. 



APPENDIX C: 

In this appendix we discuss the influence of the Pauli principle on the weight of the 
first gain satellite in the framework of time-dependent perturbation theory. This provides 
additional insight in the meaning of the result to order g A (Eq.([47|)) which we obtained from 
the exact solution Eq.(^l|) and Eqs. (f4lD and ([4^). In order to shorten the calculation and 
to simplify the argument we use the coupling function g^ l \q) which vanishes for negative q 
in this appendix and assume the phonons to be initially in their groundstate |0) ph . As we 
assume initial statistical operators corresponding to homogeneous systems, it is sufficient to 
describe the time evolution of an initial state |k ) a for N fermions and a product state for 
distinguishable particles. In the following we treat fermions and note the simplifications for 
distinguishable particles in the end. We calculate (n k ) t as 

(n k )(t) = (Mt)\Mt)} (CI) 

where \4>kif)) = n k U{t) |k ) a |0) ph with U{t) the time evolution operator in the interaction 
representation defined in Eq.(H). As we are interested in the gain region we assume that k 
is larger than the largest initial momentum k® max . For the coupling function g^'(q) the first 
order perturbation theory for U(t) gives no contribution, i.e. the leading term is of second 
order in g^ 1 ' 

/o 

-r)n k Y,l g (t)plp q \k o } a \0) ph . (C2) 
u 7 g>0 

Here we have used Eq. (f|) and the fact that the term involving p q p q > does not contribute 
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\Mt)) {2) =-"*jf dt " H e P (t')H ep (t")\k ) a \0) ph 



for k > k Q max . The function j q (t) defined in Eq.([Uj) results from the double time integration. 
Using Eq.([3|) or the representation of p q in second quantization one obtains 



h k p\ Pq |k ) Q 



i<j 

+ 



( k®, . . . , k® + q, . . . , kj q, . . . ,k 
. . . , k® — q, . . . , kj + q, . . . , ^yv) a ) 



(C3) 



where the initial momenta not written out are unshifted. Note that the states on the rhs 
can vanish if one (or both) of the shifted momenta coincide with the unshifted ones. This is 
the "trivial" manifestation of the Pauli blocking. A more subtle effect of the antisymmetry 
of the states occurs, if we perform the g-summation in Eq.(C2). This yields 

'2n- 



(2) 



L 



[7fc-fc?(^) ki, . . . , k, . . . , + kj — k, . . . , k^ 
+7fc-fc° {t) k^, . . . , + — k, . . . , k, . . . , 



(C4) 



For fermions the second term on the rhs equals the first one apart from a minus sign, 
while for distinguishable particles (without the subscript "a") they are orthogonal. Using 
Eq. (PH) we obtain for the fermionic occupation numbers in the gain regime 

-2tt\ 2 



(n k ) (2) (t) 



L 



E' 7fc-*9(f)-7*-fc<?(*) +o(r 
i<j 



(C5) 



where the prime indicates that the pair only contributes if + k® — k is different from 
all the other initial momenta kf for I different from i and j. For k — k max larger than the 
range Ak = k max — k min the prime on the sum can be omitted. For distinguishable particles 
one obtains from Eq. 

•2tt^ 2 



<»*>2 (*) 



L 



i<j 



E 7fc-fc°(*) + lk-k°(t) 



(C6) 



which also follows directly from the exact solution Eq.(p2|) summed over all particles and 
the approximation for Gq(x, t) presented in Eq. (|3^) . 

For a more quantitative discussion of the difference between the gain peaks from Eqs.( |C5| ) 
and QC5| ) it is necessary to discuss the momentum dependence of 7 9 (t) for different r = uot. 
The real part of 7 9 (t) which is equal to \f q (t)\ 2 has a peak at q B , is symmetric to this 
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momentum and apart from weak oscillations falls off like l/(q — Qb) 2 ■ The imaginary part 
of J q (t) is antisymmetric with respect to qs, with peaks near q% and falls off like l/(q — 
The width of the peaks in the real- and imaginary part of j q (t) is inversely proportional to 
t and this leads to the cancellation in Eq. (P5|) . 
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FIGURES 

FIG. 1. Evolution of the weight of the initially occupied momentum state p k o (r) in time r = uj$t 
with the coupling function g^ and the coupling a = 27rlCP 3 . As in all the following figures the 
calculations were performed with phonons initially at zero temperature and the coupling if not 
indicated differently. 

FIG. 2. Momentum distribution function for systems with different coupling g^ and g^ and 
different initial distributions at time u>ot = 50. The coupling strength is always a = 27rl0 -3 and 
the particle density n p = 0.001. For the first two systems with an initially sharp momentum 
distribution, i.e. a = the continuous part of the distribution n con t is shown with the weight of the 
initial momentum state being 0.7295 and 0.7280 for the first , and second system respectively. The 
difference between the two curves is only visible in the energy oscillations for positive momenta. 
For the third system with an initial Gaussian distribution the distribution function n (see (|A6|) ) is 
shown. 

FIG. 3. Momentum distribution function for the thermodynamic limit h con t with different 
particle densities and for a finite system with 9 particles and qs = nsAq with ng = 40 (Aq is 
the smallest momentum unit). The coupling strength is a = 2-7rl0 -2 . The particle density of the 
first two systems is the same (h p = 0.2/(27r)). The weight of the initially excited state in the 
thermodynamic limit systems is 0.040 for the first and 0.135 for the last one. For the finite system 
the zero momentum component has been left out, having the value 0.042. 

FIG. 4. Momentum distribution function for three finite systems with coupling strength 
a = 27rl0 -3 at time uj^t = 50. For the two systems with six fermions and six distinguishable 
particles we use ng = 42 and for the two fermion system ns = 14 in order to give the same 
dimensionless particle density. The occupancies of the initially occupied adjacent states at zero 
momentum have been rescaled by 1/5 to allow for a better resolution. 

FIG. 5. The six particle systems are the same as in Fig|4| for uJot = 50 but with the coupling 
a = 2vrl0- 2 . 
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FIG. 6. The systems and the time are the same as in Fig.|H| but with a = 2-7rlO -3 and different 
initial momenta. Here, three adjacent momentum states at zero momentum and three shifted by 
— qs are initially occupied. The occupancies of these states have been again rescaled by 1/5. 
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